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Abstract 

We investigate the effect of anharmonicity on the WKB approximation in a 
double well potential. By incorporating the anharmonic perturbation into the 
WKB energy splitting formula we show that the WKB approximation can be 
greatly improved in the region over which the tunneling is appreciable. We 
also observe that the usual WKB results can be obtained from our formalism 
as a limiting case in which the two potential minima are far apart. 
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It is well known that quantum tunneling leads to a splitting of degenerate energy levels 
in a symmetrical two-well potential. There are three approaches to the calculation of this 
energy splitting: the WKB approximation, the instanton method, and numerical calculation. 
From the comparison of the results from the WKB and instanton methods with those of 
numerical calculations, it was shown that the instanton method is better than the WKB 
approximation |y, because the WKB method is generally believed to have inherent errors 
associated with the connection formula ||. The modified-barrier [[| and modified- well 0] 
formalisms have been proposed for the improvement of the WKB approximation. Recently 
the authors in Ref. || have shown that a careful account of the phase changes in connection 
formula improves the accuracy of the WKB wave function. 

In this letter we propose another formalism whereby the energy splitting within the 
WKB approximation becomes consistent with the instanton result. Unlike many of the 
WKB formalisms, the present work incorporates the anharmonicity into the WKB formalism, 
which gives a more realistic model, and hence more improved energy splitting result. In other 
words, the incorporation of anharmonicity results in a level shift due to the perturbation in 
each well. 

Consider a particle of mass m in a one-dimensional symmetrical two-well potential 

2 

V{x) = r ^ r {x-af{x + af, (1) 

where u is the angular frequency in each well when the two wells are far apart, and ±a are 
the positions of the two potential minima. For a tunneling to occur the separation between 

2 2 

the two minima should be large enough so that the height of the barrier mu 8 a is higher than 
the lowest energy level in each well. In the limit a — > oo, the potential is divided into two 
independent harmonic oscillator potentials in which the lowest energies are the same and 
given by E = ^hu. When these two harmonic oscillator potentials approach each other, 
they become an anharmonic potential, so that the lowest energies are no longer Eq because 
of the anharmonic perturbation. 

To evaluate the lowest perturbation energies we expand V(x) around each minima ±a. 
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Since the potential is symmetric, we consider one of either positions. For the minimum at 
x = a we have 



mu 2 , , 9 [ x — a 3(x — a) 
- 1 + + - 



V(x) = -7— (x - a)' 



a a 2 



(2) 



Following a standard perturbation theory it is straightforward to show that the perturbation 
energy to second order correction is 

E = E Q [1 + e(rj)} , (3) 

where e{rj) is defined as 

e(7/) = || (25 - 18977 2 ), 
and we have introduced a dimensionless parameter 



mtoa? 



which is small for large a. We see that the first term in Eq.(3) corresponds to the lowest 
energy of the unperturbed harmonic potential and e{rj) is the correction term which was 
ignored in the previous studies. In the following, we demonstrate that this correction term 
plays an important role in the improvement of the WKB approximation. 
Using Eq.(3) we write the WKB level splitting formula as || 

2^ a 

AE WKB = —e~ s , (4) 



where 



S=~ I" \/2m(V(x) - E)dx, 

Yl J —a 



T = [ ^=dx, (5) 

Ja ^E-V(X) 

and ±a, ±7 are the four classical turning points (Fig. 1) corresponding to the perturbed 
energy E. a and 7 can be expressed in terms of e(rj), repectively, as 



a = adl-2 V Jl + e{r]), 7 = adl + 2 V Jl + e{rj). (6) 



Since we are interested in the region with large values of a, the elliptic integrals in Eq.(5) 
can be performed asymptotically for small rj. Keeping only the dominant terms in r/, we 
obtain the WKB energy splitting 



AE- 



WKB 
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Tiuj^—e 3^ 
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(7) 



where 5(rj) is defined as 



5{rj) 



V 71 + <V) 



exp 



£ M _ e (,) l n I ™ 1 + £(,,) 



(8) 



Comparing this with the instanton result ]7| 



AE, 



e ^ 



we find that 



WKB 



AE, 



few- 



(9) 



In the limit that the two potential minima are completely separated, which implies rj —>■ 0, 
we see from Eq.(8) that 5(rf) — > 1. In this regime the Eq.(9) reduces to 



AE^ 



(o) 



(10) 



where AE W ] KB is the WKB energy splitting obtained without considering the anharmonicity 
effect. The ratios in Eqs. (9) and (10) as a function of dimensionless parameter 7] are shown in 
Fig. 2. Note that, in the range of an appreciable tunneling probability (e.g., 0.1 < rj < 0.15), 
the WKB result obtained from the present formalism is arbitrarily close to the instanton 
result. 

A few comments are addressed in the following. Eq.(10) agrees well with the results of 
Refs.[l,2] in which the difference between the WKB approximation and the instanton method 
is claimed to be attributed to the errors introduced by the WKB connection formula. We 
note here that they |I|J^] obtained the energy splitting in the limit a — ► oo, which corresponds 
to zero tunneling probability and is not allowed in the calculation of the energy splitting due 



to tunneling. In order for the tunneling probability not to vanish, thus, the two potential 
wells should not be far apart. In this case, the coupled potential wells are simulated as 
an anharmonic potential more realistically than two idealized harmonic potentials which 
was assumed in their calculations. Our result shown in Eq.(7), which includes the effect 
of anharmonicity, is based on the formalism that includes the region where the tunneling 
probability is appreciable. As we can see in Table I, within the range of the occurence of 
a considerable tunneling, significant improvement in the usual WKB approximation can be 
achieved by the incorporation of the anharmonicity. 

In summary, we suggest a more realistic formalism with anharmonicity included than the 
previous ones only with ideal harmonicity. While the results from Refs.[l,2] are valid only in 
the limiting case of a — > oo, our approach is applicable to the broader range of the separation 
between two wells, over which the tunneling amplitude is conspicuous. Moreover, in this 
region, our formalism greatly improves the usual WKB methods, and the WKB energy 
splitting obtained from this formalism is shown to be in good agreement with that from 
the instanton method. Taking the limit a — > oo of our result, we found that our expression 
reduces to the previous one as shown in Eq.(lO) 0. Whereas it was only conjectured in 
Ref. [[J that the modification factor may come from the connection formula of the 
WKB approximation, we, here, obtained that factor through the analytical method with 
anharmonicity incorporated. 
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FIGURES 

FIG. 1. One dimensional anharmonic double well potential. E is the perturbation energy to 
second order and ±a, ±7 are the classical turning points corresponding to E. 

FIG. 2. The ratios in (9) and (10) as a function of 77 are plotted. The region of large rj has 
been excluded in this plot because the two equations (9) and (10) are asymptotic expressions for 
small 77. In the limit 77 — ► (that is, a — > 00) the two plots exactly agree with each other. 
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TABLES 



TABLE I. A comparison between AEwkb and AEi n within the range of 0.1 < rj < 0.15. 



AE WKB 

_2_ AE,„ 

0.1 0.98104 

0.121 0.99870 

0.122513 1.00000 

0.123 1.00042 

0.125 1.00214 

0.127 1.00386 

0.13 1.00644 

0.15 1.02349 
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